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It is presented a thorough analysis of scalar perturbations in the background of Gauss-Bonnet, 
Gauss-Bonnet-de Sitter and Gauss-Bonnet-anti-de Sitter black hole spacetimes. The perturbations 
are considered both in frequency and time domain. The dependence of the scalar field evolution 
on the values of the cosmological constant A and the Gauss-Bonnet coupling a is investigated. For 
Gauss-Bonnet and Gauss-Bonnet-de Sitter black holes, at asymptotically late times either power- 
law or exponential tails dominate, while for Gauss-Bonnet-anti-de Sitter black hole, the quasinormal 
modes govern the scalar field decay at all times. The power-law tails at asymptotically late times 
for odd-dimensional Gauss-Bonnet black holes does not depend on a, even though the black hole 
metric contains a as a new parameter. The corrections to quasinormal spectrum due to Gauss- 
Bonnet coupling is not small and should not be neglected. For the limit of near extremal value of 
the (positive) cosmological constant and pure de Sitter and anti-de Sitter modes in Gauss-Bonnet 
gravity we have found analytical expressions. 



PACS numbers: 04.30.Nk,04.50.-|-h 



I. INTRODUCTION 

Black holes in more than four spacetime dimensions 
are of considerable interest recently due to the two main 
reasons: they naturally appear in string theory, and in 
extra dimensional brane- world scenarios 1]. According 
to some of these scenarios it is possible that the small 
higher dimensional black holes can be produced in parti- 
cles collisions in Large Hadron Collider. At the same time 
quantum gravity may show itself already at TeV-energy 
scale. Yet, the effects of quantum gravity then may be 
observed as corrections to classical General Relativity. 

String theory predicts quantum corrections to classical 
General Relativity, and the Gauss-Bonnet terms is the 
first and dominating correction among the others. Sev- 
eral higher other theories of gravity sustain black hole 
solutions. The solution for neutral black hole in Gauss- 
Bonnet gravity was obtained by Boulware and Deser 
and Wheeler 0] . More generally. Lovelock gravity has 
been studied and shown to possess black hole solutions 
with interesting thcrmodynamical properties 0, Q . 

Thus, the problem of black hole production in trans- 
planckian particle collisions has attracted considerable 
interest recently in the context of large extra dimensions 
scenarios of TeV-scale gravity. It was observed that the 
classical spacetime has large curvature along the trans- 
verse collision plane, as signaled by the curvature invari- 
ant {RijkiK^-''^^ and thereby quantum gravity effects, and 
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higher curvature corrections to the Einstein gravity, can- 
not be ignored Q. At the same time we know that after 
formation of such a black hole its evolution has three 
stages: first, it looses its "hairs" coming into Kerr-like 
phase, then looses angular momentum transforming to 
Schwarzschild-like black hole, and finally exerts strong 
Hawking eva,£oration what results in loosing mass (see 
for instance 8] and references therein). The stage when 
black hole perturbations decay, transforming perturbed 
black hole into unperturbed one, is governed by quasi- 
normal modes and is the aim of our present research. 

In this paper we consider quasinormal perturbations 
of Gauss-Bonnet black holes including a non-vanishing 
cosmological constant. Quasinormal modes are a very 
useful tool to uncover properties of the intrinsic geome- 
try, since the modes characterizes well the geometry and 
does not depend on further extrinsic properties, inde- 
pendent of the geometry itself . They have been used 
successfully in a large class of astrophysical questions, 
from black holes to stars. In addition, it has been argued 
that the Gauss-Bonnet gravity in asymptotically anti-de 
Sitter (AdS) spacetimes may be analyzed through anti- 
de Sitter/conformal field theory (AdS/CFT) correspon- 
dence within next-to- leading order In this case the 
quasinormal modes of the large Gauss-Bonnet- AdS black 
hole could find a holographic interpretation in conformal 
field theory, as is the cases for the AdS black hole in 
Einstein gravity [Tll |. 

In |l2j the quasinormal modes for a charged asymptot- 
ically fiat black hole in Gauss-Bonnet gravity were found 
with the help of the WKB approach The tensor- 

type gravitational perturbations for Gauss-Bonnet black 
hole has been considered recently in 

To obtain the quasinormal modes we use numerical 
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analysis as well as a semi-analytical WKB-type treat- 
ment. Such an approach is based on the fact that the 
wave equation is similar to a Schrodinger equation, and 
depending on the kind of potential, it makes sense to 
borrow the methods used in quantum mechanics in or- 
der to define a semi-classical approximation. This vein 
has been followed and an approximation for the quasinor- 
mal frequencies has been obtained to a high WKB order 
[T3| . In addition to the frequency domain, we analyze 
the evolution of scalar perturbations in the time domain 
and find good agreement between the results found by 
the two approaches. 

We have observed that at asymptotically late time, 
power-law tails do not depend on the Gauss-Bonnet 
coupling and are the same as for the d-dimensional 
Schwarzschild black hole, when d is odd. For several 
simpler particular cases, namely, for pure de Sitter and 
anti-de Sitter space-time (without black hole) , and for ex- 
tremal Gauss-Bonnet-de Sitter black hole we have found 
exact analytical formulas for (quasi)normal modes. The 
QNMs for Gauss-Bonnet black holes, with coupling a ^ 1 
predicted by string theory, is seemingly different from 
those of Schwarzschild black hole. Therefore the GB- 
corrections to the QN spectrum should not be ignored, 
when considering Tev-scale of quantum gravity scenarios. 
All found modes are damping, what implies stability of 
Gauss-Bonnet black holes against scalar field perturba- 
tions. 

The paper is organized as follows: Sec. II represents 
the preliminaries of the Gauss-Bonnet- (A)dS metric and 
its scalar perturbations. Sec. Ill is devoted to the meth- 
ods used in the paper, namely the WKB method (in 
the frequency domain), and the characteristic integra- 
tion method (in the time domain). Sec. IV discuss the 
quasinormal behavior of the Gauss-Bonnet (GB), Gauss- 
Bonnet-de Sitter (GBdS) and Gauss-Bonnet-anti-de Sit- 
ter (GBAdS) black holes. In Sec. V we discuss the future 
perspective and some unsolved questions in this field. 



II. GAUSS-BONNET BLACK HOLE 
SOLUTIONS 



The Einstein-Gauss-Bonnet action in 
d-dimensional spacetime model has the form 



the 
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giRabcdR""'" - 'iRcdR'" + i?" - 2A) (1) 



where R and A are the d-dimensional Ricci scalar and 
the cosmological constant, respectively. The parameter 
a represents the (positive) Gauss-Bonnet coupling con- 
stant, which is related to the Regge slope parameter or 
string scale. 



The Gauss-Bonnet Lagrangian Cgb is given by 

Cgs = i?" - 'iR^uR'"' + Rt^.paR'"''"' ■ (2) 

One should note that in four dimensions the Gauss- 
Bonnet term ((2J) is a total divergency, and yields upon 
integration a topological invariant, namely the genus of 
the hypersurface defining the Gauss-Bonnet action (but 
even in four dimensions there are interest in the GB cor- 
rection, as seen in for example). 

A metric obtained as a solution of the field equations 
is given by 



ds^ 



-h{r)dt^ + h{r)-^dr^ + r'^dVLl_^ 



where the function h(r) is given by the expression 



h{r) ^ 1 + 



2a 




{d-l){d-2) 



(3) 



(4) 



The constant fi is proportional to the black hole mass 
and dVi^_2 is the line element of the {d — 2)-dimensional 
unit sphere. The constants a and at are connected by 
the relation: 



a = 167rGd(d- 3)(d- 4)a/ 



(5) 



We set up a scalar field $ on such a background obey- 
ing the Klein- Gordon equation 
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(6) 



In order to separate the wave function in terms 
of eigenpotential we first separate the variables as 
<^{t,r,{9i\) ^ Rf{t,r)Yf,n{{6i\)/r. As usual we obtain 
a simple equation for Ri(t,r), which is given by the ex- 
pression 



^ d^Rtiu,v) 
dudv 



+ V{r{u,v))Rt{u,v) = 



(7) 



where u = t — ri,, v = t + ri, and the tortoise coordinate 
is defined by the relation 



dr 



1 



h{r) 



(8) 



The variables u and v are the light cone coordinates cor- 
responding to the time and tortoise coordinate. The ef- 
fective potential for the scalar field in ^ is 



V{r) = h{r) 

^d~2 dh{r] 
2r dr 



{d-2){d-4) 



4^2 

£i£ + d- 



hir) 



3) 



(9) 



The effective potential is positive definite potential 
barrier for any / for GB black hole, and, for / > for 
GBdS black hole (For I = GbdS case, the negative 
pitch appears). For GBAdS case the potential diverges 
at infinity. 
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III. NUMERICAL AND SEMI-ANALYTICAL 
APPROACHES 

A. Characteristic integration 

In [T^ a simple but very efficient way of dealing with 
two-dimensional d'Alembertians has been set up. Along 
the general lines of the pioneering work light-cone 
variables have been introduced, leading to Q. 

In the characteristic initial value problem, initial data 
are specified on the two null surfaces u — uq and v — vq. 
The basic aspects of the field decay are independent of 
the initial conditions (as confirmed by simulations), so 
we use Gaussian initial conditions. 

Since we do not have analytic solutions to the time- 
dependent wave equation with the effective potentials 
introduced, one approach is to discretize the equation 
(|7|l. and then implement a finite differencing scheme to 
solve it numerically. One possible discretization, used for 
example in pL lia l2ll| . is 



Re{N) = RiiW) + R,{E) - RfXS) 



(10) 



where we have used the definitions for the points: N = 
(w A, w A), W = (m A, w), E = {u,v + A) and 
S = {u, v). Another possible scheme is 
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Re{N) = Ri{E) + ReiW) ~ Re{S) 



[ViS)ReiS) + V{E)RfXE) + V{W)Ri{W)] 



(11) 



Although the second discretization pil) is more time con- 
suming than l|10() , it was observed in 22J that Ijlll) is more 
stable for fields in asymptotically AdS geometries. With 
the use of expression H1U|I or Hll|) . the basic algorithm 
will cover the region of interest in the u — v plane, using 
the value of the field at three points in order to calculate 
it at a forth one. After the integration is completed, the 
values of Ri in the regions of interest are extracted. 



B. WKB analysis 

Considering the Laplace transformation of the Eq. {Tj) 
(in terms of t and r*), one gets the ordinary differential 
equation 

f^-[s^ + V{r{r.))]Mr.)^0 . (12) 

One finds that there is a discrete set of possible values of 
s such that the function ipilr-i,) satisfies both boundary 
conditions: 



lim ipi 

r* — *-t-oo 



(14) 



By making the formal replacement s — ico, we have the 
usual quasinormal mode boundary conditions. The fre- 
quencies Lu (or s) are the quasinormal frequencies. 

The semi-analytic approach used in this work is 
a very efficient algorithm to calculate the quasinormal 
frequencies, which have been applied in a variety of situ- 
ations p^ . 

Under the choice of the positive sign of the real part of 
uj, QNMs of Gauss-Bonnet and Gauss-Bonnet-de Sitter 
black holes satisfy the following boundary conditions 



Tpiri,) ~ C± exp(±ia;r,t) as r+ ±oo 



(15) 



lim tp£ 



r* — ^ — oo 



(13) 



corresponding to purely in-going waves at the event hori- 
zon and purely out-going waves at null infinity (or cosmo- 
logical horizon, if A > 0). For the Gauss-Bonnet-anti-de 
Sitter geometries, the effective potential is divergent at 
spatial infinity (which corresponds to a finite value of r^, 
here taken as 0). In the present work, we assume Dirich- 
let boundary conditions, setting ipi(ri, = 0) = 0. 

To find the quasinormal modes of the black hole whose 
effective potential has the form of a potential barrier (GB 
and GBdS black holes) one can use a high order WKB 
approach, finding 

^ ^ I - L2- Ls- Li~ L5- Le = , (16) 

where Vq is the height and Vq is the second derivative 
with respect to the tortoise coordinate of the potential 
at the maximum. L2, L4, and Lq are presented in 
01 . Thus wc are able to use this formula for finding the 
quasinormal modes of Gauss-Bonnet and Gauss-Bonnet 
de Sitter black holes. Yet, for Gauss-Bonnet anti-de Sit- 
ter it cannot be applied as the corresponding potential is 
divergent at spatial infinity. 

Accuracy of WKB approach may be bad for some cases 
of higher dimensional black holes. We think that it is 
mainly not because of second small peak in higher di- 
mensional case [2^ , [3^ : the WKB inaccuracy is limited 
by the case i = n or £ < n. To judge about accuracy 
of WKB method one has to compare the WKB results 
with results obtained by an accurate Frobenius proce- 
dure. This was done for a d-dimensional Schwarzschild 
black hole in a paper |3^ . where it was shown that for 
low overtones (i > n) the difference between 6th order 
WKB and Frobenius method results is less then one per- 
cent. We believe this signifies the relialability of WKB 
formulas for £ > n modes, even for higher dimensional 
black holes. After all, for Z > modes, and for scalar 
field perturbations considered in this paper, there is no 
negative pitch in the potential. 
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Figure 1: Field decay in the Gauss-Bonnet black holes, for 
d = 5 and d = 7. It is observed a quasinormal mode domi- 
nated region. Asymptotically, the field decays as a power-law 
tail (dashed lines). The parameters in this graph are a = I, 
/i = f .0 and i^O. 




Figure 2: Power-law tails in the Gauss-Bonnet black holes. 
The estimated power- law are Re oc t~^-°''^, Re oc t~^-°^^ and 
Re oc for a = 0, Of = 5 and a = 10 respectively. The 

predicted power for a = is —5. The parameters in this 
graph are d = 7, = 1.0 and £ = 0. 



IV. EVOLUTION OF PERTURBATIONS: TIME 
AND FREQUENCY DOMAIN 

In this section we shall discuss the quasinormal and 
late-time behavior for scalar field perturbations in the ex- 
terior of Gauss-Bonnet black holes, generally, with a null, 
positive and negative A-term, and therefore one has to 
consider the correlation of the scalar field evolution with 
"global" parameters: GB-coupling a, A-term, spacetime 
dimensionality, and "local" parameters such as black hole 
mass fi and multipole number £. 



A. Gauss-Bonnet black holes 

As seen in a previous work ^^l) the WKB method al- 
lows very accurate calculations of the quasinormal modes 
associated with the field evolution. A complementary 
analysis can be performed within the time-dependent pic- 
ture. For this purpose, we use here a characteristic initial 
value algorithm. 



The scenario presented by the WKB calculations is 
consistent with the results obtained with time evolution 
approach. From the wave- functions calculated with the 
characteristic integration routine, it is observed that, af- 
ter an initial transient phase, the decay is dominated by 
the quasinormal mode ringing. It is possible to estimate 
with high precision the oscillatory and exponential decay 
parameters using a non-linear fitting based in a anal- 
ysis. We emphasize that the numerical concordance is 
excellent, as seen in Tables HI - II VI The results, as com- 
pared between WKB approximation and characteristic 
integration agree to an accuracy within a few per-cents 
for Z > n case. This small difference must exist, because 
we compare the data for the fundamental overtones in 
frequency domain with time domain data where the con- 
tribution from all overtones is taken into consideration. 
Unfortunately the WKB accuracy ior I — n — is not 
satisfactory what results in large difference between fre- 
quency and time domain data for that case. 

Strictly speaking, the WKB technique we used here 
converges only asymptotically. Practically, WKB for- 
mula shows good convergence within several few orders 
after eikonal approximation. Yet, the worse convergence 
of the WKB method takes place when we deal with the 
intermediate values of a « 1. That is why, in this regime, 
the agreement between the WKB and the characteristic 
integration results is the worst. 

The imaginary part of the frequency does not show too 
much dependence on a, yet slightly decrease when a is 
increasing. On the other hand, the real part increases 
with a, though not significantly either. This might be 
showing that a quasinormal mode is much more an ef- 
fect connected with the local geometry containing the 
black hole rather than with the global effect of the ge- 
ometry, namely, the effect of the existence of an event 
horizon matters much more than a detailed dependence 
on the parameter a. Yet, large enough values of a cer- 
tainly affect the quasinormal spectrum: the QNMs are 
proportional to a in the regime of large a As a 

approaches zero, the QNMs go to those of ordinary d- 
dimensional Schwarzschild black hole described by the 
Tangherlini metric. 

However, it does not mean that GB corrections are 
negligible. On the contrary, according to string theory 
the Gauss-Bonnet coupling a should be around 1. Let 
us compare the results, for instance, for £ = 2, n = 
QNMs for Schwarzschild and a = 1 GB black holes for 
d — 6: for Schwarzschild we have uj = 1.5965 — 0.3967i 
(6th order WKB), for a = 1 (6th order WKB) we get uj = 
1.69654 — 0.31929i (Note that for this case convergence is 
good and the 3th order WKB value is not much different 
1.67624- 0.323698i). Thus the effect of GB coupling is 
about 6.3% in real and more then 20% in the imaginary 
part here. For larger values of a it is certainly larger. 
We have the same order of difference for other values of 
n and £. 

In the time domain the signal has three stages: the ini- 
tial pulse dependent on the source of perturbations, the 



5 



Table I: Values for the quasinormal frequencies for the funda- 
mental mode in the Gauss-Bonnet geometry, obtained from 
sixth order WKB method and directly from characteristic 
data, for d = 5, ^ = 0, 1 and several values of a. 

d = 5 WKB Characteristic Integration 

i a Re(aJo) -lm{u)o) Re(ci;o) -Im(tt;o) 



Table IV: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet geometry, obtained from 
sixth order WKB method and directly from characteristic 
data, for d = 8, ^ = 0, 1 and several values of a. 

d = 8 WKB Characteristic Integration 

£ a Re((j;o) -lm{ijJo) Re(LJo) -lm{u)o) 






0.1 


0.389935 


0.256159 


0.379 


0.282 





0.2 


0.396034 


0.250548 


0.383 


0.272 





0.5 


0.429741 


0.208293 


0.391 


0.246 


1 


0.1 


0.720423 


0.255506 


0.7234 


0.250 


1 


0.2 


0.723177 


0.252684 


0.7284 


0.245 


1 


0.5 


0.739205 


0.236885 


0.7443 


0.228 



Table II: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet geometry, obtained from 
sixth order WKB method and directly from characteristic 
data, for d = 6, ^ = 0, 1 and several values of a. 



d = 6 WKB 

£ a Re(ti;o) -lm(tJo) 



Characteristic Integration 
Re(a;o) -Im(a;o) 



0.1 0.735854 0.402416 0.7109 



0.2 

0.5 

5 

10 

20 



0.748053 
0.83530 
0.906661 
1.513624 
2.961675 



0.391049 
0.304837 
0.144820 
0.456935 
0.927128 



0.7144 
0.7225 
0.9526 
1.5182 
2.878 



0.1 1.139007 0.415034 1.153 



1 0.2 

1 0.5 
1 5 

1 10 

1 20 



1.136193 
1.158635 
1.790103 
3.260415 
6.437139 



0.415706 
0.391339 
0.258382 
0.498426 
0.991374 



1.159 
1.176 
1.791 
3.253 



0.412 
0.402 
0.375 
0.226 
0.423 
0.952 

0.395 
0.386 
0.363 
0.253 
0.504 



Table III: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet geometry, obtained from 
sixth order WKB method and directly from characteristic 
data, for d = 7, ^ = 0, 1 and several values of a. 



d = 7 WKB 

£ a Re(ti;o) -Im(ti;o) 



Characteristic Integration 
Re(ci;o) -lm(ajo) 



0.1 
0.2 
0.5 

5 
10 
20 

0.1 
0.2 
0.5 

5 
10 
20 



1.11738 
1.13699 
1.29469 
1.39823 
1.48053 
2.00368 

1.54573 
1.53194 
1.56277 
2.01379 
2.47824 
3.39094 



0.546056 
0.529543 
0.395111 
0.574472 
0.385616 
0.56458 

0.577608 
0.58701 
0.554551 
0.308316 
0.243737 
0.59452 



1.092 
1.092 
1.092 
1.275 
1.515 
2.001 

1.587 
1.530 
1.609 
1.982 
2.475 
3.387 



0.532 
0.520 
0.493 
0.351 
0.358 
0.567 

0.527 
0.517 
0.489 
0.337 
0.423 
0.597 






0.1 


1.51702 


0.694245 


1.461 


0.676 





0.2 


1.54463 


0.673566 


1.463 


0.658 





0.5 


1.7854 


0.488086 


1.469 


0.616 





5 


1.47941 


0.868859 


1.647 


0.420 





10 


1.800 


0.410252 


1.838 


0.421 





20 


2.15426 


0.556954 


2.154 


0.544 


1 


0.1 


1.93407 


0.750046 


2.021 


0.652 


1 


0.2 


1.90391 


0.773947 


2.024 


0.637 


1 


0.5 


1.94463 


0.734591 


2.035 


0.602 


1 


5 


2.44511 


0.455166 


2.348 


0.434 


1 


10 


2.6644 


0.463966 


2.670 


0.458 


1 


20 


3.2276 


0.581024 


3.209 


0.590 



quasinormal ringing dominating period, and the power- 
law tail (see Fig^. The bigger GB-coupling is, the 
larger the quasinormal dominated region, i.e. at later 
times the tails start dominating. As can be seen from 
Figl^l the power-law tails do not show dependence on 
the Gauss-Bonnet coupling a and are the same as for the 
d-dimensional Schwarzschild black hole in Einstein gen- 
eral relativity, when d is odd. That is, the fields always 
shows a power-law falloff: for odd d > 3 the field behaves 
as 



(17) 



at late times, where £ is the multipole number. This 
behavior is entirely due to d being odd and does not 
depend on the presence of a black hole [3| ■ It is known, 
that in Einstein gravity, for even d > 4, the field decays 
as ^ 



Ri OC t' 



-{2i+3d-S) 



(18) 



and for the latter case there is no contribution from the 
flat background. This power-law tail is entirely due to 
the presence of the black hole j2j|. At the same time, 
the Gauss-Bonnet black hole metric H3I4|) goes to pure 
Minkowskian metric when the black hole mass equals 
zero, i.e. in a space-time without a black hole. In other 
words, empty space-time in Gauss-Bonnet gravity "does 
not see" the a. That is why we do not observe the a- 
dependence of tails in odd space-time dimensions. 

Thus, if the Gauss-Bonnet term changes late-time be- 
havior, it must show itself only for even-dimensional 
space-time. In the numerical procedure developed with 
the characteristic integration scheme, no tails (power-law 
or otherwise) were observed in even dimensions Gauss- 
Bonnet spherical black holes. Yet, it should be pointed 
out that the integration of the scalar field equation in 
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the GB background is a much more demanding numer- 
ical problem than the same integration with the usual 
Einstein coupling. In the latter case there are auxiliary 
analytical results, such as explicit expression for the tor- 
toise coordinate function. Therefore, the GB codes are 
less precise and more time consuming, and eventual tails 
could be hidden. The possible absence of tails with even 
d deserves further consideration. 



B. Gauss-Bonnet-de Sitter black holes 

For the GBdS black holes the quasinormal ringing 
stage becomes correlated with a new parameter: a pos- 
itive A-term. When the A-term is growing, both real 
oscillation frequency and the damping rate are decreas- 
ing. Yet, real part of uj is more sensitive to the changes 
of A-term. 

Qualitatively this resembles the quasinormal oscilla- 
tions of d-dimensional Schwarzschild-de Sitter black hole 
[2^. In the limit of extremal value of the A-term, i.e. 
when the cosmological horizon (rc) is very close to the 
event horizon {rj-), it is possible to generalize the for- 
mulas found in [2g for four dimensional black hole and 
in |23| for (i-dimensional case. Namely, the quasinormal 
frequencies for the near extreme Gauss-Bonnet asymp- 
totically de Sitter black holes are given by 



li{e + d-3)rc-r+ 1 
n = 0,1,2,... , 



(19) 



where (a function of a in this generalized context) is 
the surface gravity at the event horizon. The above for- 
mula is well confirmed numerically: Fig|3| shows a com- 
parison of the values obtained by direct numerical calcu- 
lation and from Eq. H19(l . Thus the quasinormal modes 
are proportional to the surface gravity at least for 
lower overtones. It should be pointed that for the usual 
Schwarzschild-de Sitter black holes, numerical and ana- 
lytical investigations [2^ suggest that the high overtone 
behavior does not obey the formula ((T^ . 

When the mass parameter ^ is set to zero, we have 
the case of pure de Sitter spacetime in the Gauss-Bonnet 
gravity. The metric function 10} then reduces to the fol- 
lowing form: 



^^^ = ^+2^-2^1 



8aA 



(d-l)(d-2) 



(20) 



Repeating the analysis of |2g|,|30|, we come to the conclu- 
sion that quasinormal modes exist only in odd spacetime 
dimensions and are given by the formula: 





J , \ , \ , I , I I 

20 40 60 80 100 
1(1 + 5) 



Figure 3: Graph of Re(aj)^ as a function of £{£ + d — 3), in 
the near-extreme limit positive A limit. The bullets are the 
values calculated from the time-evolution profiles, and the 
dashed lines are values obtained from the expression 1191 . The 
parameters in this graph are k+ — 0.01, a — 1, jj. — 1.0, and 
the differences between the analytical and numerical results 
are under 2%. 



Note that pure Gauss-Bonnet-de Sitter quasinormal 
modes are purely imaginary, which corresponds to ex- 
ponential decaying without oscillations. 

It is well-known that the late-time tails of black holes 
in asymptotically de Sitter space-time for zero multipole 
and for higher multipoles are qualitatively different. For 
the zero multipole field {£ = 0), the time domain picture 
is the following: after a transient part, a quasinormal 
mode dominated region is best observed. Following the 
quasinormal mode dominated region, a late-time decay 
region settles. In this latter phase, the wave-functions 
decay asymptotically to a constant value, as has been the 
case in the Schwarzschild dc Sitter black hole which was 
studied before 18, 19, 31). This is illustrated in Figgl 

For first and higher multipoles {£ > 0) at late times we 
observe exponential tails in vicinity of GBdS black hole. 
This is also an expected result, since in Einstein gravity 
the exponential tails are observed as well in usual de Sit- 
ter black holes [HllllIJl. This is illustrated in FigEl 
The dependence of the quasinormal modes on A-term and 
Gauss-Bonnet coupling can be learnt from Tables X-XII 
for different space-time dimensionality. 

The numerical simulations developed for Gauss- 
Bonnet-de Sitter black hole indicate that the massless 
scalar perturbation in this geometry behave asymptoti- 
cally as 



i?f « exp c9""'^K^) as t->oo 



d> 4 and a > 



(22) 



where Hc is the surface gravity at the cosmological hori- 
zon and c'""'' is an adjustment parameter for the cor- 
rection. The expression H22|l shows that, although the 
exponential tail in the GBdS background is dependent 
on the Gauss-Bonnet coupling a (since function 
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Table V: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet-de Sitter geometry, ob- 
tained from third and sixth order WKB method, for d = 5, 
£ — and several values of a and A. 



d = 


5 


WKB (3th order) 


WKB (6th order) 


a 


A 


Re{ujo) 


-Im(L<;o) 




-lm(wo) 


0.1 


1/8 


0.30485 


0.278407 


0.334981 


0.26202 


0.1 


1/2 


0.14922 


0.217542 


0.152444 


0.2164 


0.1 


2/3 


0.0677006 


0.14792 


0.00666957 


0.150988 


1 


1/5 


0.301869 


0.218752 


0.378865 


0.159611 


1 


1 


0.112613 


0.148749 


0.11619 


0.145465 


1 


7/5 


0.0240051 


0.0683652 


0.0242716 


0.0709775 



Figure 4: Field decay in the Gauss-Bonnet-de Sitter black 
holes, with £ = 0. A quasinormal mode dominated region 
is observed (above), and asymptotically the field decays to 
a constant (below). The fundamental mode, calculated with 
the WKB method and directly from the characteristic data 
are 0.8356 - 0.2935i and 0.8011 - 0.2608i. The parameters in 
this graph are d = 6, a = 3.0, A = 0.1 and n = 1.0. 




Figure 5: Field decay in the Gauss-Bonnet-de Sitter black 
holes, for d = 5,6, 7, 8. It is observed a quasinormal mode 
dominated region. Asymptotically, the field decays as an ex- 
ponential tail. The parameters in this graph are a — 0.1, 
H = 1.0 and £ ^ 1. 



of a), the form of the dependence is identical to the null 
a case. Eg generalizes the analogous expression 

found in "Sri for the usual Schwarzschild black holes. 



C. Gauss-Bonnet-anti-de Sitter black holes 

The quasinormal and late-time behavior of black holes 
in anti-de Sitter spacetime is significantly different from 
those in asymptotically de Sitter or flat spacetimes. The 
key difference is stipulated by the effective potential be- 
havior, which is divergent at spacial infinity. Thus the 
anti-de Sitter space acts as an effective confining box. 
Therefore the Dirichlet boundary conditions are natu- 
ral. These boundary conditions are required also by 
AdS/CFT correspondence for scalar field perturbations 
[Tl|. Yet, for higher spin perturbations the true bound- 
ary conditions may be different 20] . 

In the usual Schwarzschild- anti-de Sitter black holes, 
the quasinormal modes govern the decay at all times and 



Table VI: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet-de Sitter geometry, ob- 
tained from third and sixth order WKB method, for d = 6, 
£ — and several values of a and A. 



d = 


6 


WKB (3th order) 


WKB (6th order) 


Q 


A 


Re(tJo) 


-lm(tjo) 


Re((x)o) 


-Im(ci;o) 


0.1 


1/4 


0.580472 


0.428367 


0.652079 


0.411337 


0.1 


1 


0.379931 


0.392027 


0.401925 


0.384688 


0.1 


2 


0.0653851 


0.168544 


0.0636328 


0.173794 


1 


1/2 


0.582293 


0.358225 


0.736989 


0.235326 


1 


2 


0.318697 


0.307084 


0.331268 


0.282461 


1 


4 


0.00390515 


0.0487306 


0.0574555 


0.0529852 


10 


100 


1.5033 


0.574511 


1.5741 


0.599309 



Table VII: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet-de Sitter geometry, ob- 
tained from third and sixth order WKB method, for d = 7, 
£ — and several values of a and A. 



d = 


7 


WKB (3th order) 


WKB (6th order) 


a 


A 


Re(LJo) 


-Im(c<;o) 


Re(a;o) 


-lm(wo) 


0.1 


1 


0.769109 


0.556647 


0.852395 


0.551612 


0.1 


1 


0.561644 


0.517639 


0.598484 


0.511069 


0.1 


1 


0.00963007 


0.222564 


0.0939457 


0.228191 


1 


1 


0.868246 


0.48557 


1.12603 


0.2999562 


1 


4 


0.461426 


0.419016 


0.483098 


0.381285 


1 


7 


0.113334 


0.215329 


0.111663 


0.218394 



Table VIII: Values for the quasinormal frequencies for the 
fundamental mode in the Gauss-Bonnet-de Sitter geometry, 
obtained from third and sixth order WKB method, for d = 8, 
£ — and several values of a and A. 



d = 8 




WKB (3th order) 


WKB (6th order) 


a 


A 


Re(ajo) -lm{ujo) 


Re((.Jo) -lm{ujo) 


0.1 


1 


1.1449 0.669364 


1.29739 0.695376 


0.1 


4 


0.644892 0.598472 


0.682218 0.594199 



8 



Table IX: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet-de Sitter geometry, ob- 
tained from third and sixth order WKB method and directly 
from characteristic data, for d — 5, £ — 1 and several values 
of a and A. 

d = 5 WKB (6th order) Characteristic Integration 



a 


A 


Re ((.Jo) 


-lm(ajo) 


Re{ujo) 


-Im(a;o) 


0.1 


1/8 


0.64317 


0.241723 


0.6451 


0.2382 


0.1 


1/2 


0.379265 


0.165101 


0.3856 


0.1564 


0.1 


2/3 


0.231024 


0.103394 


0.2247 


0.1051 


1 


1/5 


0.698655 


0.181854 


0.6926 


0.1877 


1 


1 


0.365998 


0.111668 


0.3704 


0.1133 


1 


7/5 


0.153506 


0.0465747 


0.1518 


-0.042631 



Table XI: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet-de Sitter geometry, ob- 
tained from third and sixth order WKB method and directly 
from characteristic data, for d = 7, ^ = 1, a = 0.1 and several 
values of A. 

d = 7 WKB (6th order) Characteristic Integration 



a 


A 


Re(a;o) 


-lm(tJo) 


Re(ci;o) 


-Im(a;o) 


0.1 


1 


1.28621 


0.519229 


1.304 


0.4910 


0.1 


2 


1.00125 


0.440876 


1.008 


0.4306 


0.1 


4 


0.303405 


0.153651 


0.3091 


0.1416 


1 


1 


1.48546 


0.411346 


1.467 


0.3623 


1 


4 


0.899964 


0.338928 


0.8944 


0.3224 


1 


7 


0.361602 


0.152214 


0.3638 


0.1398 



Table X: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet-de Sitter geometry, ob- 
tained from third and sixth order WKB method and directly 
from characteristic data, for d — 6, £ — 1 and several values 



of a 


and A. 








d = 


6 


WKB (6th order) 


Characteristic Integration 


a 


A 


Re(tJo) 


-Im(a;o) 


Re((^o) 


-lm(ajo) 


0.1 


1/4 


1.04698 


0.39859 


1.0574 


0.3828 


0.1 


1 


0.74644 


0.323514 


0.7479 


0.3183 


0.1 


2 


0.238066 


0.114496 


0.2387 


0.1212 


1 


1/2 


1.09756 


0.301244 


1.053 


0.2948 


1 


2 


0.6962 


0.244326 


0.6962 


0.2307 


1 


4 


0.0847573 


0.0315107 


0.08965 


0.03323 


10 


100 


3.21004 


0.52559 


3.1940 


0.5360 


10 


2000 






2.1711 


0.4900 


10 


5000 






0.6228 


0.2018 



thereby no power-law or exponential tails appear |21ll22| . 
We observe a similar behavior the scalar field perturba- 
tions in the Gauss-Bonnet-anti-de Sitter black holes. 

It is not possible to use WKB method to find the quasi- 
normal modes in Gauss-Bonnet- AdS case because the ef- 
fective potential is not a potential barrier anymore. The 
Horowitz-Hubeny method [TJ is not applicable either, 
because the Taylor expansion of the effective potential 
has infinite number of terms. That is why we were lim- 
ited only by time domain analysis, which is free of the 
above problems. From FigEl we see that, indeed, the 
quasinormal modes are dominating even at sufficiently 
late times. We also have learnt from FigElthat the quasi- 
normal mode dominated region grows, as the multipole 
index £ grows. 

As is known from Einstein action the radius of 

the AdS black hole goes to zero, the quasinormal modes 
of the black hole approach its pure anti-de Sitter values 
[3^ . Repeating the calculations of j3^] , we find the exact 



Table XII: Values for the quasinormal frequencies for the fun- 
damental mode in the Gauss-Bonnet-de Sitter geometry, ob- 
tained from third and sixth order WKB method and directly 
from characteristic data, for d = 8, ^ = 1, q = 0.1 and several 
values of A. 



d = 


8 


WKB (6th order) 


Characteristic Integration 


a 


A 


Re(ajo) -lm(LJo) 


Re((^o) -lm{ujo) 


0.1 


1 


1.75036 0.697491 


1.803 0.6250 


0.1 


4 


1.10954 0.508301 


1.103 0.4964 



expression for the normal modes in GB gravity: 

1/2 




(2n- 



(23) 



The pure GB-AdS modes, unlike GB-dS modes, exist in 
any any spacetime dimension. 



V. CONCLUSIONS 

We have considered here frequency and time domain 
description of evolution of scalar field perturbations in 



Table XIII: Values for the quasinormal frequencies for the 
fundamental mode in the Gauss-Bonnet-anti-de Sitter geom- 
etry, estimated from the characteristic data, for d = 5, ^ = 0, 
fi = 1.0, A = —0.1 and several values of a. 



Re(cjo) 



-Im((x)o) 



0.1 
0.1 
0.5 
1.0 
1.5 



0.4923 
0.4920 
0.4904 
0.4885 
0.4866 



-0.01585 
-0.01593 
-0.01634 
-0.01702 
-0.01766 
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Figure 6: Field decay in the Gauss-Bonnet anti-de Sitter black 
holes, for several values of the Gauss-Bonnet coupling. It is 
observed a quasinormal mode dominated region. Asymptoti- 
cally, the field decays in quasinormal modes. The parameters 
in this graph are d = 5, = 1.0 A = —0.1 and £ = 0. 



the exterior of black holes in Gauss-Bonnet theory of 
gravity, generally with a A-term. The quasinormal be- 
havior even though being corrected by a new parame- 
ter, Gauss-Bonnet coupling a, are qualitatively depen- 
dent mainly on the A-term and black hole parameters 
such as mass ^ and multipole number £. The late-time 
tails for asymptotically flat Gauss-Bonnet black holes, do 
not depend on the Gauss-Bonnet coupling in odd space- 
time dimensions, and therefore are the same as those 



for d-dimensional Schwarzschild black hole in Einstein 
gravity. Moreover, in the case of Gauss-Bonnet-de Sitter 
black holes, the late-time tails, though dependent on a, 
yet, rather trivially, i.e. only through dependence of the 
surface gravity at the cosmological radius on a. Thus, 
the Gauss-Bonnet coupling shows itself "minimally" in 
late-time behavior. The most interesting problem which 
remains unsolved is, to find late time tails for even di- 
mensions, and thereby, to know whether the power-law 
tails depend upon the Gauss-Bonnet term. At the same 
time, we have shown that corrections to the quasinormal 
frequencies due to GB-term are not negligible: they may 
reach 20% for string theory predicted values of a ~ 1. 

Even though our analysis can easily be extended to the 
massive scalar field, we were limited here by the mass- 
less case. We expect that the influence of the massive 
term upon the QNMs will be similar to that found in 
[33| . i.e. the lower overtones should be corrected by the 
field mass, infinitely high overtone asymptotic will be un- 
changed no matter the value of the massive term. Also 
we did not consider the high overtone behavior of the 
GB black holes. Generally, the high overtone asymj 
totics must be studied by totally different methods 
and deserves separate investigation. 

Acknowledgments 



This work was partially supported by Fundagdo de 
Amparo d Pesquisa do Estado de Sao Paulo (FAPESP) 
and Conselho Nacional de Desenvolvimento Cientifico e 
Tecnologico ( CNPq), Brazil. 



[1] L. RandaU and R. Sundrum, Phys. Rev. Lett. 83 3370, 
4690 (1999). 

[2] D. G. Boulware and S. Deser, Phys. Rev. Lett. 55, 2656 
(1985). 

[3] J. "Wheeler, Nucl. Phys. B268, 737 (1986). 

[4] D. Lovelock, J. Math. Phys. 12, 498 (1971). 

[5] J. Crisostomo, R. Troncoso and J. Zanelli, Phys. Rev. D 
62, 084013 (2000). 

[6] Elcio Abdalla and L. Alejandro Correa-Borbonet, Phys. 
Rev. D 65, 124011 (2002). 

[7] V. S. Rychkov, Phys. Rev. D 70, 044003 (2004). 

[8] G. Duffy, C. Harris, P. Kanti, E. "Winstanley, 
hep-th/0507274 (2005). 

[9] K. D. Kokkotas and B. G. Schmidt, Living Rev. Relativ- 
ity 2, 2 (1999). 

[10] S. Nojiri, S. D. Odintsov, J. of High Energy Phys. 0007, 
049 (2000). 

[11] G. T. Horowitz and V. E. Hubeny, Phys. Rev. D 62, 
024027 (2000); A. Nunez, A. O. Starinets, Phys. Rev. D 
67, 124013 (2003). 

[12] R. Konoplya, Phys. Rev. D 71, 024038 (2005) 
hep-th/0410057 . 

[13] B. F. Schutz and C. M. WiU, Astrophys. J. 291, L33 
(1985); S. Iyer and C. M. "WiU, Phys. Rev. D 35, 3621 
(1987); R. A. Konoplya, Phys. Rev. D 68, 024018 (2003) 



|gr-qc/030305 2l. 

[14] G. Dotti, R. J. Gleiser, gr-qc/0503117' 

[15] Rodrigo Olea, J. of High Energy Phys. 0506, 023 (2005). 

[16] C. Gundlach, R. Price and J. Pullin, Phys. Rev. D 49, 

883 (1994). 
[17] R. Price, Phys. Rev. D 5, 2419 (1972). 
[18] P. R. Brady, C. M. Chambers, "W. Krivan and P. Laguna, 

Phys. Rev. D 55, 7538 (1997). 
[19] P. R. Brady, C. M. Chambers, "W. G. Laarakkers and E. 

Poisson, Phys. Rev. D 60, 064003 (1999). 
[20] I.G. Moss and J.P. Norman, Class. Quant. Grav. 19, 2323 

(2002) 

[21] B. "Wang, C. Molina and E. Abdalla, Phys. Rev. D 63, 
084001 (2001). 

[22] Bin "Wang, Chi-Yong Lin, C. Molina, Phys. Rev. D 70, 
064025 (2004). 

[23] E. Berti, K. D. Kokkotas, Phys. Rev. D 71, 124008 
(2005); R. A. Konoplya, Phys. Rev. D 66, 084007 (2002) 
gr-qc/0207028 ; S. Fernando, Gen. Rel. Grav. 37, 585 
(2005); R. A. Konoplya, E. Abdalla, Phys. Rev. D 71, 
084015 (2005) hep-th/0503029 ; R. A. Konoplya, C. 
Molina, Phys. Rev. D 71, 124009 (2005) gr-qc/0504139 '; 
H. Nomura, T. Tamaki, Phys. Rev. D 71, 124033 (2005). 

[24] V. Cardoso, S. Yoshida, O.J.C. Bias, J. P.S. Lemos, Phys. 
Rev.D 68, 061503 (2003). 



10 



[25] R. A. Konoplya, Phys. Rev. D 68, 124017 (2003) 

hep-th/0309030 . 
[26] V. Cardoso and J. P. S. Lemos, Phys. Rev. D 67 084020 

(2003). 

[27] C. Molina, Phys.Rev. D 69, 104013 (2004). 

[28] S. Yoshida and T. Futamase, PRD 69, 064025 (2004); 
V. Cardoso, J. Natario, R. Schiappa, J. Math. Phys. 45 
4698 (2004); R. A. Konoplya and A. Zhidenko, J. of High 
Energy Phys. 037, 06 (2004) hep-th/0402080 . 

[29] E. Abdalla, K. H. C. Castello-Branco, A. Lima-Santos, 
Phys. Rev. D 66, 104018 (2002). 

[30] J. Natario, R. Schiappa, hep-th/0411267 (to be published 
in Adv. Theor. Math. Phys.). 

[31] C. Molina, D. Giugno, E. Abdalla, A. Saa, Phys.Rev. D 
69, 104013 (2004). 



[32] R. A. Konoplya, Phys. Rev. D 66, 044009 (2002) 

hep-th /0205142 . 
[33] k. A. Konoplya, Phys. Lett. B 550, 117 (2002) 

"gr-qc/0210105 ; R. A. Konoplya, A. V. Zhidenko, Phys. 

Lett. B 609, 377 (2005) gr-qc/0411059 . 
[34] H-P NoUert, Phys. Rev. D 47, 5253 (1993); V. Cardoso, 

R. Konoplya, J. P.S. Lemos, Phys. Rev. D 68, 044024 

(2003) gr-qc/0305037 ; L. Motl and A. Neitzke, Adv. 
Theor. Math. Phys. 7, 307 (2003). 

[35] E. Berti, M. Cavagha and L. Guahieri, Phys. Rev. D69, 
124011 (2004); H. Yoshino, T. Shiromizu and M. Shibata, 
gr-qc/0508063 

[36] V.Cardoso, J. Lemos, S. Yoshida, Phys. Rev. D 69, 044004 

(2004) . 



